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Abstract 



We introduce a class of randomly time-changed fast mean-reverting stochastic volatility 
(TC-FMR-SV) models. Using spectral theory and singular perturbation techniques, we derive 
an approximation for the price of any European option in the TC-FMR-SV setting. Three 
Qh , examples of random time-changes are provided and are shown to induce distinct implied 

volatility surfaces. The key features of the TC-FMR-SV framework are that (i) it is able to 
incorporate jumps into the price process of the underlying asset (ii) it allows for the leverage 
effect and (in) it can accommodate multiple factors of volatility, which operate on different 
time-scales. 



^ ■ 1 Introduction 

Stochastic volatility models have played an important role in the derivatives markets over the past 
(NT twenty years. Much of the success of stochastic volatility models is due to the fact that two of 

the earliest and most well-known models-the Heston model [T2] and the Hull- White model |23j- 
capture the most salient features of the implied volatility surface while simultaneously preserving 
the analytic tractability needed to quickly calculate the price of an option. Yet the short-comings 
of these models is well-documented in literature. For example, the Heston model misprices far in- 
and out-of-the-money European options [111136] . 

There are a number of possible explanations for why the earliest stochastic volatility models 
fail to match implied volatility levels across all strikes and maturities. One theory is that a single 
1_h 1 factor of volatility, running on a single time scale, is not sufficient for describing the dynamics of 

the volatility process. Indeed, the existence of several factors of volatility has been documented 
in literature [2]|U|HJ[TU][Tni[20 ,25,26, 28 . Such evidence has led to the development of multi-scale 
stochastic volatility models, i.e. models in which instantaneous volatility levels are controlled by 
multiple diffusions running of different time-scales [T3"l[T7ll3T)] . 

Another line of reasoning states that jumps in the underlying asset price are required in order 
to capture the true dynamics of the market. Empirical work supports this notion [S]. Hence, 
academics and practitioners have developed models that incorporate both jumps in the asset price 
as well as stochastic volatility [51151 |3~I]. 

Along these lines, Mendoza-Arriaga et al. recently introduced a unified credit-equity framework 
in which the underlying asset is modeled as a stochastically time-changed scalar diffusion [27] . This 
work is notable for a number of reasons. First, the scalar diffusion that controls the asset price 
may exhibit both local volatility (i.e. volatility that is a function of the scalar diffusion itself) and 
killing (i.e. jump to default). When the local volatility is modeled as a negative power of the scalar 
diffusion a decrease in the underlying asset price results in an increase in volatility. This feature, 
known as the leverage effect, has been empirically documented [7J. Additionally, by subjecting the 
scalar diffusion to a random time-change the authors are able to incorporate jumps in the asset 
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price as well as non-local factors of stochastic volatility. Finally, it is shown that, under relatively 
benign conditions, the framework of Mendoza-Arriaga et al. remains analytically tractable. We 
see great value in the work of Mendoza-Arriaga et al. and seek to build upon it. 

In this paper, rather than base our model upon a scalar diffusion as in |27) . we begin with the 
class of fast mean-reverting stochastic volatility (FMR-SV) models considered by Fouque et al. 
in |14j . Such models are important because they capture the empirically known-to-exist short time- 
scale of volatility [15 ,21. Additionally, FMR-SV models capture the leverage effect by negatively 
correlating the Brownian motions that drive the asset price and volatility processes. Using the 
methods outlined by Mendoza-Arriaga et al. in [27], we subject the FMR-SV class of models to 
a random time-change. For certain classes of time-changes this has the effect of adding jumps 
to the underlying asset price as well as additional factors of volatility. These additional factors 
of volatility operate on a different time-scale than the fast mean-reverting factor volatility. We 
refer to this class of models as the class of time-changed fast mean-reverting stochastic volatility 
(TC-FMR-SV) models. 

The rest of this paper proceeds as follows. In section [2] we introduce the class of TC-FMR-SV 
models. This is done in a few steps. First, in section [2~T1 we review the class of FMR-SV models 
considered in [14] . Next, in subsection 12. 2\ we explain how the FMR-SV class can be extended 
using random time-changes. Finally, in section 12.31 we review the three classes of random time- 
change. Some specific model assumptions are listed in section |2~41 

In section [3] we develop our option-pricing methodology. Again, this is done in several steps. 
First, in section [37X1 we review some important results from spectral theory, which we immediately 
apply to the European option-pricing problem in the TC-FMR-SV setting. This reduces the 
option-pricing problem to that of solving a single eigenvalue equation. In section 13.21 we find 
an approximate solution to this eigenvalue equation using techniques from singular perturbation 
theory. Then, in section 13.31 we show how to relate the approximate solution of the eigenvalue 
equation to the approximate price of a European option. The main result of our work is the 
formula we provide in Theorem[4]for the approximate price of a European option in the TC-FMR- 
SV setting. 

In section [4] we prove the accuracy of our option-pricing approximation. And in section [5] we 
provide examples of four different random time-changes and calculate the approximate price of a 
European call option in these time-change regimes. 

2 Model Framework 

In this section we introduce a class of TC-FMR-SV models. We begin by reviewing the FMR-SV 
class of models considered by Fouque et al in [T4] . 

2.1 Review of FMR-SV Models 

Under the physical measure P the FMR-SV class of models has the following dynamics 

X = x, (1) 

Y e = y, (2) 

Here, Wt and B t are Brownian motions under P with instantaneous correlation p G [—1,1]. The 
process St represents the price of a non-dividend paying asset (stock, index, etc.), which has 
expected geometric growth rate /3 > and stochastic volatility f(Y t e ) > 0. The process Y t e 
appears as an Ornstein-Uhlenbeck (OU) process with long-run mean m £ R and "vol of vol" 
v > 0. The OU process operates on time-scale e > 0, which is intended to be small so that the 



S t =cxp (pt + Xt), 
dX t 



\f{Y{)dt + f(Yf)dW u 



dYf = -(m- Y t e ) dt + 

e Ve 

d(W,B), = pdt. 
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rate of mean-reversion (1/e) of Y t e is high. It is in this sense that Y t e is fast mean-reverting. In 
fact, Y t e need not be an OU process specifically. The essential aspect of Y t e is that it be an ergodic 
process with a unique invariant distribution. The function f(y) is left unspecified, as only certain 
moments of f(y) play a role in the FMR-SV framework. Specific assumptions on the function 
f{y) and the process Y t e will be given in section \2 .41 The parameter e will play an important role 
throughout this paper. As such, we will use a superscript e to indicate dependence on this small 
time-scale parameter. 

For the purpose of option-pricing, it is necessary to move to the risk-neutral pricing measure, 
which we denote as P. Under P the FMR-SV class of models has the following dynamics 



St = cxp(rt + X t ) , 
dX t = ~f\Yf)dt + f{Yf)dW u 

<iy; = I(m-y/)-^r(y/) 



d{W,Bj = pdt. 



dt + -^dB t , 



X = x, 



Y a = y, 



(3) 



(4) 



Here, Wt and B t are Brownian motions under P with instantaneous correlation p. The Girsanov 
transformation, which relates the physical measure P to the risk-neutral measure P, is chosen such 
that the volatility-driving process Y t e acquires a market price of volatility risk T(Y t t ) and such 
that the discounted asset price (e~ rt St) is a martingale under P. Note that r > is the risk-free 
rate of interest. As was the case with f(y), the function T(y) is left unspecified, as only certain 
moments of T(y) play a role in the FMR-SV framework. Specific assumptions on T(y) will be 
given in section 12^41 



2.2 TC-FMR-SV Models 

Under the risk-neutral measure P, the TC-FMR-SV models have the following form 

S t = exp(ri + A Tt ). 

Here, (X t ,Y t e ) is as described by equations ([S]) - Q in section l2Tl The key difference between 
the TC-FMR-SV class of models and the FMR-SV class is that the dynamics of the log of the 
discounted asset price log (e - St), which would simply be given by the two-dimensional Markov 
diffusion {X t , Y^) in the FMR-SV framework, is now given by a time-changed diffusion {Xx t , Y^ ). 
Broadly speaking, a random time-change T t is an increasing process, starting from zero, which is 
independent of (Xt,Y t e ) and is right-continuous with left limits. We list specific assumptions on 
the random time-change T t in section \2. 41 



2.3 Stochastic Time-Changes 

In this paper, we will consider three classes of stochastic time-changes: Levy subordinators, abso- 
lutely continuous time-changes, and time-changes that are the composition of a Levy subordinator 
and an absolutely continuous time-change. These three classes of stochastic time-change are em- 
ployed extensively in 27 in the context of local volatility models with state-dependent killing 
rates. Drawing inspiration from [27], we will use these classes of time-change in the context of 
FMR-SV models. A review of each of these classes is presented below. In an effort to avoid 
re-inventing the wheel, our discussion will be brief, focusing mainly on those aspects necessary for 
calculating option prices. For a more detailed discussion of stochastic time-changes, we refer the 
reader to [27] . 
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2.3.1 Levy Subordinator T/ 

A Levy Subordinator is a non-decreasing Levy process with positive jumps and non-negative 
drift. Because all Levy processes have stationary and independent increments, the Laplace trans- 
form of a Levy subordinator can be expressed as 



E 



e -ATi 



= e -0(A)t 



(A el), (5) 



Ael:E 



< 00 



The function <f>(A) is known as the Levy exponent of the subordinator and is given by the 
Levy-Khintchine formula 



0(A) = 7 A+ / (l-e- As )is(ds). (6) 
Jo 

Because all Levy subordinators are of finite variation no truncation of integral © is necessary. 
The absence of a A-independent constant term in ^ means that we have excluded any killing of 
the stochastic time-change. We require that the drift 7 of the subordinator T t be non-negative 
7 > 0. 

The Levy measure v, which must satisfy 

(1 A s) v{ds) < 00, 

describes the arrival rate and distribution of jumps. Specifically, for some Borel set B G 23 (R+), 
the value v(B) gives the intensity of a Poisson process that counts the number of jumps of size 
s G B. 

For A > expectation ([5]) is always finite. However, in order to prove the accuracy of our 
pricing approximation in section |4j we will need to consider the case A < 0. To characterize the 
set I, we recall Theorem 25.17 of [33], where it is established that 



E 



<oo V( <^ / e~ As v{ds) < 



In general, I is an interval (A, 00) or [A, 00) where A < 0. 

An important sub-class of Levy subordinators are the subordinators of compound Poisson 
type. The jump component of such subordinators is described by a compound Poisson process 
with (net) jump arrival intensity a > and jump size distribution F. For such subordinators, the 
Levy measure u(ds) can be written 

v(ds) = aF(ds), 

in which case the Levy exponent, given by equation ([6]), becomes 

0(A) = 7A + a ( 1 - / e As F(ds)) . (7) 



Although it is not strictly necessary for our framework, for the sake of computational simplicity, 
we will be primarily interested in Levy subordinators for which the Levy exponent <fi(A) is known 
in closed form. 

We would like to emphasize the importance of Levy subordinators as a class of stochastic time- 
change. Because Levy subordinators T t exhibit jumps, the time-changed diffusion (X T i , K^i ) (and 
thus the asset price St) will exhibit jumps as well. To our knowledge, this is the first time that 
jumps in the asset price St have been incorporated into the FMR-SV framework. 

This concludes our brief review of Levy subordinators. For more thorough coverage, we refer 
the reader to [6]. 
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2.3.2 Absolutely Continuous Time-Change T 2 



We now consider stochastic time-changes of the absolutely continuous type. When we say T 2 is 
an absolutely continuous time-change we mean that T 2 can be written as 



T 2 = / V(Z s )ds, 
Jo 



Z = z, 



(8) 



where Z t is an infinite lifetime Markov process taking values in Mr. The function V : M — > R+ 
shall be referred to as the rate function of stochastic time-change. We are primarily interested in 
absolutely continuous time-changes T 2 for which the Laplace transform 

L(t,z,A) =E Z [e- AT * 2 ] , (9) 
iJe- A /» y H, (A G Jt) 



Jt--=\Ae 



E z 



< OO 



is known explicitly. Here, the notation E 2 [•] is used to indicate the conditional expectation 
E[-\Z = z\. We note that © is always finite for A > 0. Values of A < for which © is 
finite depend on the specific choice of T 2 and must be checked on a case-by-case basis. 

In the previous section, we showed that Levy subordinators are an important class of stochas- 
tic time-change because jumps in the subordinator induce jumps in the asset price. Absolutely 
continuous time-changes are important for a very different reason; they have the ability to change 
a one-factor stochastic volatility model into a multi-factor stochastic volatility model. To see 
this, we define (^X t ,Y^j :— ^AT T 2,Y^ 2 ^. Then, there exist P-Brownian motions W t and B t with 
correlation p such that [29 



dX t = ~-f 2 (Yf)V(Z t )dt + f(YnvnzT)dw t , 



dY e 



7^ 



r(n € 



iV2 



V(Z t )dt + ^y/V(z7)dB t , 



X a = x, 

Yo = y- 



Note that the volatility of X t is controlled by the product f(Y^)yJV{Z t ) rather than just the single 
factor f(Y t e ), which controls the volatility of X t . Note also that the multiple factors of volatility 
are operating on different time-scales; f(Y t e ) acts on a time-scale of 0(e) and y / V(Z t ) acts on 
a time-scale of 0(1). As demonstrated in [13], when compared to their one-factor counterparts, 
multi-factor stochastic volatility models in which the factors of volatility operate on different 
time-scales have the ability to vastly improve the fit to the empirically-observed implied volatility 
surface. 



2.3.3 Composite Time-change T 4 3 

Finally, we may consider composite time- changes, which are time-changes of the form 

Here, T} is a Levy subordinator and T 2 is an absolutely continuous time-change, which is inde- 
pendent of Tf. As long as the Levy exponent 0(A) of and the Laplace transform L(t, z,A) of 
T 2 are known explicitly, the Laplace transform of the composite time-change Tf can be calculated 
as well. This is accomplished by conditioning on the absolutely continuous time-change T 2 as 
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follows 



-ATf 



E, 



E„ 



-ATI, 



-0(A)T t 2 



= i(t,2,0(A)), 

Kt '■= |a G M : E^ 



(A e AC t 



(10) 



-AT/ 



< 00 



Again, we remark that [0,oo) C fCt- However, values of A < for which A g K, t depend on the 
specific choice of time-change Tf . 

The importance of composite time-changes is as follows: by combining a Levy subordinator 
with an absolutely continuous time-change we are able to incorporate jumps in the asset price St 
as well as add multiple factors of volatility to the class of FMR-SV models. The variety gained 
by combining different types of stochastic time-changes provides us with considerable modeling 
flexibility. 

Remark on Notation 

Throughout this paper we shall use the superscripts 1, 2, 3 to specify which type of random time- 
change we wish to consider. The notation will be used to denote a Levy subordinator, the 
notation T t 2 will be used to denote an absolutely continuous time-change and the notation T t 3 will 
be used to denote a composite time-change. Finally, if we do not wish to specify a particular class 
of random time-change we will omit the superscript altogether and use the notation Tt. 

2.4 Specific Model Assumptions 

We have now described the TC-FMR-SV class of models. However, we have not been specific 
about certain technical assumptions on the process Y t e , the functions f(y) and r(y) and the 
stochastic time-change T t . The purpose of this section is to list these assumptions in one place. 
Our assumptions are as follows: 

1. Under the physical measure P, the process Y" t e is ergodic and has a unique invariant distri- 
bution Fy, which is independent of e . We note that this implies the moments E[|y t e | fc ] are 
uniformly bounded in t. That is, for every k £ 1 + there exists a positive constant Ck < oo 
such that 



supE y [|F t e ' fc 



]<C k . 



2. The volatility function f(y) is a strictly positive function such that 

(a) the process (X t ,Y t e ) as the strong solution to SDE's ([TJ and ([2]) exists and is unique 
under P, 

(b) Jp(y)F Y (dy) <oo, 

(c) a solution $(y) to the Poisson equation ([29)) exists and is at most polynomially growing. 



3. The function T(y), which describes the market price of volatility risk, is such that 

(a) the process (X t ,Y t e ) as the strong solution of SDE's © and (UJ) exists and is unique 
under P, 

(b) there exists a finite constant Cr > such that |r(y)| < Cr- 



(:> 



4. The random time-change T t is a strictly increasing cadlag process, which is independent of 
(X t ,Y t e ) and satisfies Tq = and 



Cn 



T t e 



T t Cl 



< 00. 



2.5 The Martingale Condition 



Although we specified the class of TC-FMR-SV models under a supposedly risk-neutral measure P, 
we have not yet shown that the non-dividend-paying asset in these models satisfies the martingale 
condition 



E[e- rt *S t2 \F tl ] =e~ rtl ,S tl , 



(ii < t 2 ) . 



which is required in order for P to actually be risk- neutral. In fact, because (x T i , J , \X T 2 , Y^ 2 , Z t 

and (^X T 3, y^ 3 , Z^j are time- homogeneous Markov processes, as rigorously established in [27] . the 
martingale condition reduces to 



(11) 



where we have used the short-hand notation E a: .j, )Z [-] to denote the conditional expectation E[-|Xo = 
x, Yq — y, Zq — z]. We can verify equation (jlip by conditioning on the random time-change T t as 
follows 



%, v ,» [e~ rt S t ] 



exp - ~ £ f{Y*)dt + £ f{Y?)dw)j 



E 



x,y,z 



I 1 

exp \ x — — 



f(Y*)dt 



T, 



f(Y t e )dW t ) 



where we have used the fact that fa f 1 ( Y t) dt +S f( Y t) dW t j s an exponential martingale. Having 
established that the discounted stock-price process (e~ rt St) is a martingale under P in the TC- 
FMR-SV framework, we now move on to the option-pricing problem. 



3 Option Pricing 

In this section, we discuss how the approximate price of any European option can be calculated 
in the TC-FMR-SV setting. 

3.1 Spectral Representation of European Option Prices 

We begin with a brief review of some important results from spectral theory and semigroup 
operators. 

Theorem 1. Suppose C is a self-adjoint operator acting on a Hilbert space H.. Consider the 
eigenvalue equation for (-C) 

-Cip\ = Xijj\. 

We denote by E the projection-valued spectral measure of (— C) and by u{—C) the spectrum of 
{—£■)■ Then: 
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1. The operator (— C) has the following spectral representation 

-C= I XE(dX). 

M-c) 

If g is a real-valued Borel function on M., then g(C) can be defined via operational calculus 
and is given by 

g(-£)= [ g(X)E(dX). 

2. If there exists a number 7 > —00 such that A > 7 for all A £ cf(— L), then 

Qt = e- t( -- c \ (0<i<oc), 

defines a strongly continuous one-parameter semigroup and the solution to the Cauchy prob- 
lem 

{—d t + C)u = 7 u(0, x) — u (x), {uqE%), 

is given by 

u{t 7 x)=Q t u (x)= [ e~ xt E{dX)u (x). (12) 

Ja(-C) 

Proof. Items [1] and [2] are classical results from functional analysis. The proof of item [1] can be 
found in chapter 8 of [31] and the proof of item [2] is given in Chapter 13 of [32]. □ 

For convenience, we will write (|12p as 

/ e' xt E{d\)u {x) - f e- xt Mx)Vu (d\), (13) 

where fi uo (d\) = (dip\,u ). This will help to make it clear that E (dX) u (x) is a projection of 
uq(x) onto the eigenspace H\ :— {ip G H : — Lip — Xip}. In the special case when the spectrum 
<f(— C) is purely absolutely continuous with respect to the Lebesgue measure, equation p3]) can 
be written 

J C u] e-^ t ^ u] (x)doj, 

where the C u are constants chosen such that J C u ip u {%) du) = uq(x). 
Now, consider a function u £ (i, x, y), defined as 

u e {t,x,y) := E [h(X t )\X = x,Y £ = y] , 
The backward variables x, y satisfy the Kolmogorov backward equation 

(-d t + C^y) = 0, (14) 

u e (0,x,y) = h(x). (15) 

Note that the dt term carries a minus sign because t is a forward variable. We use the notation 
Cxy to indicate the infinitesimal generator of the Markov process (X t ,Y t e ), defined in ([3]) - (j4]). 
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For clarity, we write C x Y explicitly and state its domain T> (£ x Y ) 

£x,y = \({m- y) d y + v 2 d 2 yy ) + -L (pvy/2f{y)dl y - v^T^d, 



1 



1 



s.t. lim 



^x, v \g{X u Y t )] -g(x,y) 



exists for all (x, y) € 



Now, suppose we have the solution to the following eigenvalue equation 



-x,y*a 



A- 



(16) 



By, solution to the eigenvalue equation we mean that we have the full set of eigenvalues {A c } and 
corresponding eigenfunctions {^ A (x,y)} for which (|16[) holds. Then if the operator C XY were 
self-adjoint on some Hilbert space H, by Theorem [TJ the solution to Cauchy problem (TH1) - (fT5j) 
could be expressed as 



u e {t,x,y) 



,-A E t 



(17) 



where the measure pit would be determined by BC (1151) . 

It is not apparent that there exists a Hilbert space W on which C x Y is self-adjoint. As such, it 
is not clear at this point that eigenvalue equation (|16p has a solution, nor is it clear that u e {t, x, y) 
has a representation of the form (|17[) . Nevertheless, in this paper we do not endeavor to solve 
the full the eigenvalue equation (|16l) . Rather, we shall use techniques from singular perturbation 
theory to find an approximate solution to (fl6|) . We will show that the operator associated with the 
lowest order solution to (fl6|) is in fact self-adjoint on some Hilbert space. As a result, u e (t,x,y) 
can be approximated by a function of the form (|17|) . For the moment, however, it will ease our 
calculations if we assume that u e (t,x,y) can be written as ([TT|). which we shall refer to as the 
spectral representation of u e (t, x, y). 

Supposing u e (t, x, y) can be approximated by a function of the form p7|) . we would like to use 
this knowledge to find a spectral representation for the price of a European option in the TC- 
FMR-SV framework. To this end, we consider a European option with payoff h(St) at maturity 
date t < oo. Using risk- neutral pricing, and the Markov property of (-Xr t , Yf, , Zt), we may write 
the price of a European option P e {t, x, y, z) as 



P e (t,x,y,z) =e- rt E XiVi 
Conditioning on the random time-change T t we find 



[hi 



rt+X T , 



)} 



P f (t,x,y,z) 



h(e rt+x ^) 



Ti 



[u e (T t ,x,y;t)} 



(18) 



Note that t is just a parameter here-not a variable of u e (T,x,y;t). Now, we use ([T7|) to replace 
u e (T,x,y;t) with its spectral representation. We have 



P e (t,x,y,z) = e 



e- AeT ^ A (x,y)fi h (dA<;t) 



= e 



E- 



-A'T t 



y A (x,y)n h (dA e ;t), 



(19) 



where passing the expectation through the integral is allowed by Fubini's theorem. We have used 
the notation /x/j (eiA e ;£) to remind us that Hh(dA e ]t) depends on the paramter t through the BC 
u e (0,x,y;t) — h(e rt+x ). Assuming it exists, we refer to (|19l) as the spectral representation of the 
option price P e (t,x,y, z). Note that the expectation E 2 [e _A£T *] is given explicitly by either ([5]), 
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(p?J) or (jlOp . depending on the type of random time-change. Hence, in order to fully specify the 
price of the option P e {t, x, y, z), what remains is to solve eigenvalue equation (1161) and determine 
the measure \ih (rfA e ; t) from the BC 

U c(0,x,y;t) = h(e rf + x )- (20) 
3.2 Asymptotic Analysis of the Eigenvalue Equation 

For general f(y) and T(y) there is no analytic solution to the eigenvalue equation = C x Y ^% + 
A 6 ^^. However, we note that C x Y can be conveniently decomposed in powers of ^Je as 

C XY = lc(-2)+l £("!)+ £(0), 

e ye 
C { -V = {m-y)d v + v*dl v , 

= puV2f(y)d 2 xy - iV2r(y)fy, (21) 

£(0) = _l /2(y) ^ + I/2( y )^. 

This decomposition suggests a singular perturbative approach. To this end, we expand ^> A and 
A e in powers of i/e. We have 

^ = #i° ) +V^i 1) +e*i 2) + --- ) (22) 
A e =A(°> + ^A« + eA( 2 ) + .... 

Expanding in powers of ^fe (rather than some other power of e) is a natural choice given the 
form of C x Y . The validity of this expansion will be justified in section 01 when we establish the 
accuracy of our pricing approximation. 

We now insert the expansions for ^ A (x,y) and A e into eigenvalue equation (TlBl and collect 
terms of like-powers of \/e. The 0(e~ 1 ) and ©(e" 1 / 2 ) equations are 

0{e- x ): = £(- 2 %i 0) , 

Oie- 1 ' 2 ): = £H^m + £(-l)*(°). 

Noting that all terms in C^ 2 ^ and £( _1 ) take derivatives with respect to y, we may (and do) choose 
solutions of the form ^ A = ^ A (x) and ^ A — ^ A \x) (i.e. functions of x only). Continuing the 
asymptotic analysis, the order 0(e°) and ©(e 1 / 2 ) equations are 

O(e ) : = £(- 2 %f + + A< >) *<°> , (23) 

0( e i/2) : = £(- 2 )*f +£(- 1 )*f + (£(°) + A< >)*£ ) +A< 1 >*£ ) , (24) 

where we have used JCS- -1 ^^ (x) — in ([231 . Equations (f!23"]) and (pHf . respectively, are Poisson 

(2) (3) 

equations for ^ A (x, y) and "i 1 A ( x j 2/) m the variable ?/ of the form 

= £ Y ^ + g (25) 

where C Y = C^~ 2 ^ is the infinitesimal generator of under the physical measure P. We wish to 
consider only those solutions ^(y) of (|25|) that exhibit at most polynomial growth as y — > ±oo. 
With this restriction, a necessary condition for the solvability of (P25l) is 

(g) ■= I g(y)F Y (d y ) = 0. (26) 
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We remind the reader that Fy is the invariant distribution of Y t e under the physical measure P. 
Equation (1261) is referred to as the centering condition. Please refer to Appendix [X] for a treatise 
on the Poisson equation and the centering condition. Throughout this paper, the notation (•) will 
always indicate averaging with respect to the invariant distribution Fy. In equations (|23p and 
(|24| the centering conditions become 



0= ((£(°))+A(°))*i 0) , (27) 
= (tf-^f ) + + A<°>) + A«*i 0) - (28) 

Eigenvalue equation (|27p can be solved explicitly, as the operator (C^) is given by 

(£ (0) ) = y(t-^). - 2 :=(/ 2 )- 

However, in order to solve equation (|28|) . we need an expression for (^C^ 1 ^^ (x, -)^. To this 
end, we note from (f23|) 

£(-2)^(2) = _ ^(0) + A (0)^ $ (0) = _ ^(0) _ ^(0)^ $ (0) = _1 (/2 _ ct2) _ ^ ^(0) 

Now, introducing $(y) as a solution to the following Poisson equation 

£(" 2 )$ = / 2 - a 2 , (29) 

we may express ^ A (x,y) as 

*i 2) (z, y) = -\m {dL d x ) *f\ x ) + c(x), 

where C(x) is some function which is independent of y. Hence, using (|21[) we find that (tf-VV^fa- 
is given by 

(ct-VvPix, •)) = ^puV2f(.)dl y - uV2T(.)d y ) (C - d x ) < 0) (x) + C(xfj \ 

= A^^(x), 

where 

= V 3 {dl xx - C) + V 2 (d 2 xx ~d x ), V 2 = JL (Td y $) , V 3 = (/a s $) . (30) 
Thus, from (1281) we have 



= .4« *f + ((£W) + A<°>) + A^tff • (31) 

Given a solution to (|27|) . one can use (|3 1 1) to find expressions for ^^'(x) and A^. 

This concludes our asymptotic analysis of eigenvalue equation (|16p . Before we present an 
explicit solution to (|2T|) and (|3"Tj) . we recall the following result from Sturm-Liouville theory. 

Theorem 2. The eigenf unctions ^P'(x) of equation (|27l) /orra a complete orthonormal basis in 
the Hilbert space % := L 2 (R, s(x)dx) where 



s(x) dx = e x dx, (u,v) s — / u(x) i>(a;) s(x) dec 



TTie notation u(x) indicates the complex conjugate ofu(x). 



1 We note that H29H satisfies the centering condition and &(y) exists by assumption I2cl of section f2. 41 
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Proof. The proof is by showing that (C^) of equation (|27p is self-adjoint in L 2 (R : s(x)dx). This 
is a standard result of Sturm-Liouville theory. Details can be found in any number of texts on 
differential equations [1,3, 22, 351157] . 

Theorem 3. The order O (e°) eigen] "unctions ^^\x) and eigenvalues a£ are given by 



A.™ 



a 

T 



,(iu!+l/2)x 



where T/ie order O (e 1 ^ 2 ) corrections ^J^\x) and are 

A« = -V, 



1\ 3 /. l\ 2 \ „ /V. l x 2 



- iw + 



-y 2 



— \ioj + 



□ 

(32) 
(33) 

(34) 
(35) 



where Vi and V3 are defined in (I30|) . M^e nofe i/iai A^ ' > A^] n 



Proof. A direct substitution shows that (|32|) . (I33|) . (|34|) and (j35|) satisfy equations (|27|) and 
(f3"T|) . One can easily verify that the O (e°) eigenfunctions | \f'£ ^ (x) j form a complete basis in 
L 2 (R, s(x)dx) and satisfy the orthogonality condition 



(36) 
□ 



3.3 Option Prices 



We have found explicit expressions for the approximate eigenvalues A^ w Ai?' + v / eAL 1 ' 1 and 
approximate eigenfunctions \l>^(x, y) ~ (x) + y'e'^w (z)- We now use these expressions to 
specify the approximate price P e (t,x,y, z) w p(°>(t, x, z) + \feP^ x \t,x, z) of an option. The 
following Theorem serves as the main result of our work. 



Theorem 4. TTie approximate price of an option is given by 
P £ {t,x,y,z) « P (0) (t,x,z) + ViP (1) (t,x,z), 



-A'f'T, 



P(°\t,x,z) =e~ rt J C^(t)E 

pW(t,x,z)=e- rt [ C^(t)E z [(-A«T t 

where the coefficients cffl (t) are given by 

^) = (*L 0) (-),Me rt + -) 
for a composite time-change Tf we have 

E z [< 
E z \(-A^T^ 



-M 0) T t 



1), 



= L(f,z,0(A( J o) ) 
= d a L (t, z, 4> 

= 0(aL o) ), 

= 9 a 0(A«a) 



(37) 
(38) 

(39) 

(40) 
(41) 



a=A£ '/AL 1) 
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The corresponding expressions for a Levy subordinator T t x and an absolutely continuous time- 
change T t 2 can be recovered by setting L(t,z,cp) = e~^ t and 0(A) = A respectively. The O (e°) 

eigen) 'unctions ^^\x) and the approximate eigenvalues AJ, rs Ao, + s/e a£ are given in Theorem 

m 

Proof. Consider the spectral representation of u e (T, x,y;t) given by dnj). Recall that we use the 
notation fXh {dA^t) to remind us that it e (T, x, y; t) has a BC u e (0, x, y; t) = h(e rt+x ) that takes t 
as a parameter. We expand /i^ (dA^t) and e~ A ^ T in powers of ^fe 

li h (dAl ■ t) = C<°> (t) duj + C W (t) + . . . , (42) 
e -ALT = e ~A<°»T + ^ (_ A «t) e - A ^ 0,T + . . . . (43) 

Note that we have expanded the measure fih (<^A^;t) in terms of a density C^(t) as the spectrum 
of the (e°) eigenvalue problem is absolutely continuous with respect to the Lebesgue measure 
duj. Inserting expansions (f22l) . (|42)) and (|43)) into (JTTJ) and collecting terms of like-powers of i/e 
yields 

O(e ) : uW{T,x;t) = J (t) e^^H^, 

o(e^) : u«(2>;t) = J (c« (i) e^^) + C^{t) (-A«t) e-^tfgV)) 

where we have dropped the C^\t) e~ A( ° >T ^/^\x) term because ^/^\x) — 0. Expressions for 
cL Q) (t) and C$\t) can be obtained from the BC's u(°'(0,i;i) = h{e rt+x ) and u^(0,x;t) = 0. 
We have 

O (e°) : m (0) (0, X] t) = h(e rt+x ) = J (t)*£ 0) (x)du, 

O (e 1 ' 2 ) : u« (0, s; t) = = f tf« (t)*£°) (z)du;. 



Hence 



O (e°) : (*< 0) (-),Me rt + ')) s - _/ C<°)(i) (*C°) , *£ )) ^ da; = Cf)(i), 

O ( £ 1/2 ) : = / <#>(*) (flf, *«>)*,, = c p(t), 



where we have used (|36j) . 

We have now obtained an explicit expression for u e (T, x, y; t) (T, x;t) + -y/e (T, x; t). 

In order to find an expression for the approximate price of an option P e (t, x, y, z) w P*- ' (£, x, z) + 
i/cP^(t,x,z) we simply insert our expansion for u e (T t ,x,y;t) into (|18j) . which yields (|37|) and 
(I3"8l) . Expressions (|4T)f and (j4"Tj) are given for a composite time-change T t 3 and can be obtained by 
expanding L (t, z, 0(A^)) in powers of y/e. □ 



Remark. We note that the existence and finiteness of expectations (|40[) and (|4Tj) is guaranteed for 
all w G R by assumption 0] of section HOI and by the fact that Vw £ 1 we have — a£ ' < — A^< n < 

cf. 

Corollary 5. The function y/eP^(t,x,z) is linear in the group parameters 

Proof. First, we note that V% and Vs do not appear in C<i (t), A^ or ^^\x). Next, from ([35]) we 
see that pW(t, x, z) is linear in A£ , which, from (|35|) , is a linear function of Vz and V3. Corollary 
[5] follows immediately. □ 
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Corollary [S] relates to a very important feature of the TC-FMR-SV pricing methodology. Con- 
sider first the FMR-SV framework (no time-change). To describe a particular model within the 
FMR-SV class, one would have to specify an ergodic diffusion Yf, a market price of volatility 
risk r(y) and a volatility function f(y). For the purposes of illustration, we chose to specify Y t € 
as an OU process. This choice led us to introduce five unobservable parameters (to, e, v, p, y). 
Note however, that neither the value of these parameters nor the precise form of T(y) and f(y) 
are required in order to calculate the approximate price of an option (the approximate price of 
an option in the FMR-SV framework is given by setting T t = t in the TC-FMR-SV framework). 
Rather, to O(yje), the approximate price of an option can be expressed in terms of (<r 2 , V^, V£) 
as well as the observable parameters (t,x,r). 

As mentioned previously, the particular choice of V t e as an OU is not central to our analysis. 
We could have simply written Y t € under the physical measure P as 

where a{y) and (3(y) are such that the assumptions of section I2T41 are satisfied. In this case, the 
group parameters would have become 

vi = (Pfdyfy , v 2 £ = ^~ (prd v $) . 

The key point is that, when the volatility-driving process Y t e is fast mean- reverting and satisfies 
the conditions of section HOl the details of the process are unimportant. In terms of option-pricing, 
to 0(y/e), all that matters are the values of (a 2 , V 2 e , V3 ). This is true regardless of the particular 
choice of Y t e . 

In the TC-FMR-SV framework the situation remains the same - to calculate the approximate 
price of an option, precise knowledge of the volatility-driving process Y t e is not required. However, 
the particular choice of random time-change T t does affect the approximate price p(°'(t, x, z) + 
y/e P^\t, x, z) of an option. Thus, when calibrating a particular model within the TC-FMR-SV 
class to fit market data (be the data quoted option prices or implied volatilities), the unobservable 
parameters that must be extracted are (a 2 , V 2 £ , V£) as well as the parameters of the random time- 
change T t . We will show in section [5] that different time-changes induce distinct implied volatility 
surfaces. Whether the introduction of time-change parameters is justified by the modeling flexi- 
bility the random time-change provides is a topic left for future research. 

4 Accuracy of the Approximation P e w + y/e 

In the previous section, we gave a derivation of the approximate price of a European option 
pe _ p (o) + yip(i) 

using singular perturbative arguments. The purpose of this section is to 
establish the accuracy of this approximation. In addition to the assumptions listed in section 12. 41 
we shall need one additional assumption for our accuracy proof. 

• The payoff function h(e rt+x ) and all derivatives taken with respect to x are smooth and 
bounded. 

Obviously, the most common options - calls and puts - do not fit this assumption. To prove 
the accuracy of our pricing approximation for calls and puts would require regularizing the option 
payoff as was done for the class of FMR-SV models in [16]. The regularization procedure is beyond 
the scope of this paper. As such, we limit our analysis to options with smooth and bounded payoffs. 
Before stating our main accuracy result we need the following Lemma. 

Lemma 6. Suppose J(y) is at most polynomially growing. Then, there exists a constant C such 
that 

& y [J(Y t e )} < Ce tc ?. 
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Proof. First, we define M tl the exponential martingale used in Girsanov's theorem to transform 



the measure on Y t e from P to 



M t := cxp (- J T{Y^)dB s - 1 jf P 



We note 



[Af 2 ] = E, 
<E„ 
= e * c rE 



exp ( - jf (2r(y/)) dB a -\j^ (2r(y/)) 2 d s + jf (r 2 (r/)) ds 
cxp (- jf* (2T(i7)) ds s - i jf (2r(y/)) 2 d s + jf c 2 d s 



exp - / 2r(57)dB a 



2 Jo 



(2r(r/)) 2 d s 



where we have used assumption I3bl of section R2T41 to bound T 2 (Y t e ) by C 2 . Hence 



[\Y?\ k ] = E y [\n\ k M t ] < (l„ [|F/| 2fc ] E, [A/ 2 ] ) 1/2 < (C 2k e tc < 



1/2 



The first inequality is an application of Cauchy-Schwarz. The second inequality follows from the 
above bound as well as assumption [T] of section I2T41 Since J(y) is bounded by a polynomial, this 
proves lemma |6l □ 

From here, we shall proceed as follows. First, we shall establish the accuracy of the approxi- 
mation u e (T,x,y;t) s» u^(T, x; t) + y/t vP-> (T, x; t). Then, we show how this result can be related 
to the accuracy of the approximate option price P e (t,x,y,z) sa p(°>(t,x, z) + y/eP^(t,x,z). 

Theorem 7. For fixed (T, x, y, t) there exists a constant C such that for any e < 1 the solution 
u e (T,x,y;t) to PDE {H} with BC ® satisfies 

u e - (u^ + yfeu^l <eCTe TC l 

Proof. First, we define a remainder term i? e (T, x, y; t) 

i? e := u e - U°) + yfeu^ W 2) + ev W 3 >) . 

Next, we see that 

(-dr + C e x y) ^ - -C^u^ - -L ( C(-*>uV> + tf-VuW) 
e y/e \ J 

- (c<--*>u<V + C^u^ + (-dr + r^) uW) 

- V~e (C^u^ + + (-d T + £<°>) U W) 

- e ^-^^ + (-a T + £ (0) ) « (2) + \/i (-Sr + £ (0) ) u (3) 



-dr + £x y) R 



d T + £% Y )R e = -eF', 



F e := 



(tf-Vu® + (-&r + £ (0) ) w (2) + V~e (-dr + £ (0) ) w (3) ) 
i? £ (0,x,y;i) = eG e (x,y;t), 
G e (x,y;t) :=- (u< 2 >(0,s,y; t) + V£«< 3 > (0,s,y;t)) . 



(44) 
(45) 
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Now, from the Feynman-Kac formula we note that R e (T,x,y;t), which is the solution to PDE 
with BC P5)l . has the following stochastic representation: 



R e (T,x,y;t) = eE x , 



G e (X T , Yt\ t)+ F e (s, X s , Y s e ; t) ds 



As established in [14], from the boundedness of the payoff function h(e rt+x ), and from the as- 
sumptions of section T2.41 one can deduce that F e (s,x,y;t) and G e (x,y;t) are bounded in x and 
at most polynomially growing in y. Hence, by lemma [6j there exists a constant C\ such that 



|i? e (T,x,y;t)| <eCiTe i0 r. 



Therefore, 



< |i? e |+e U ( 2 > + ^ 3 ) 

< eCiTe^ +eC 2 

< eCTe TC r, 

for some constants C2 and C. This establishes Theorem [7] 
Now we state our main accuracy result. 



□ 



Theorem 8. For fixed (t, x, y, z) there exists a constant C such that for any e < 1 the price of a 
European option P e {t,x,y, z) given by (|18[) . satisfies 



P e - 



< eC. 



Proof. 



P e (t, x, y, z) - (t, x, 2) + V~eP {1) (i, a, z)) 

u £ (T t , x, y, z; t) - (u (0 )(T t , x, z; i) + u (1) (T t , x, 2; t)) 



e~ rt E 



x.y.z 



E, 



E, 



; (T t , x, y, 2; t) - («(°) (T t , x, 2; i) + (T t , x, 2; t)) 



<E T 



< E T 



u e (T t ,x,y,z;t) - [u {0 > (T t ,x, z;t) + y/eu { ' (T t ,x, z;t) 



<edC T 
= eC, 



eCiT t e TtC ? 



(by Theorem [7]) 

(by assumption 2] of section 



for some constants C\ and C. This establishes Theorem [5] 



□ 



5 Call Option Examples 

In this section we provide examples of how to calculate the price of a European call option in 
four different time-change regimes. These examples demonstrate both the flexibility and analytic 
tractability of the TC-FMR-SV framework. 
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5.1 FMR-SV 



The first regime we consider is that of no random time-change (i.e. T t = t). This choice for T t 
reduces the TC-FMR-SV framework to that of pure FMR-SV. To calculate the approximate price 
of a European call option p(°\t,x) + ^/eP^(t,x) we use equations (l37|) and (|38[) of Theorem |4j 
Since T t is not random in the present scenario, expectations (|4"0")) and (|4"Tj) reduce to 



_ A (0) Tt l _ -A<°>i 



The O (e°) eigenfunctions *L 0) (a;) are given in (j3"2")l . Expressions for A^ -* and Ai, 1 -* can be found 
in (j3"3")l and ([55)1 respectively. Hence, what remains in order to calculate the approximate price of 
a call option P(°)(i, x) + \/e P^\t, x) is an expression for Ci 0) (i). 

For a European call with strike price K = e k and time of maturity t < oo, the option payoff 
h(S t ) is given by 

h(S t ) = (S t -e k ) + . 



Using equation (|39|) we calculate 



= f -^e(~ iu +^ 2 > (e rt+x - e k ) + e- x dx (46) 
Jk v27r 

e fe(l/2-iw) u iuJ _ 2) - e rt+A(l/2-iw) (4^ + 2 ) 

= === -. (47) 

VMl+4w 2 ) v ' 

Note that integral (|46[) will not converge for purely real values of uj. However if we move uj into 
the complex plane uj = u r + ioJi and we fix the imaginary part of u) such that Ui < (—1/2), then 
integral (|4*5]) will converge. Upon doing this, when calculating option prices using (|57)l and (|55]l. 
we must remember to hold the imaginary part u>i < (—1/2) fixed and integrate with respect to 
the real part of u (i.e. set dw = dtu r ). 

Figured] demonstrates the implied volatility surface induced by the TC-FMR-SV framework in 
the T t = t regime. We plot implied volatilities I versus log- moneyness-to- maturity ratio (LMMR) . 
We remind the reader that / and LMMR are defined by 

P BS (t, K, I) = (P (0) + V^P {1) ) («, K), 
LMMR = log (K/S ) It, 

where P BS {t, K, I) is the Black-Scholes price of a call option with strike price K , time to maturity 
t and volatility /. The notation (P^ - 1 + y/ePW) (t,K) is used here to indicate the approximate 
price of call option as calculated in the TC-FMR-SV framework with strike price K and time to 
maturity t. The parameters used in figure [T] are 

r = 0.00, a = 0.34, V 2 e = 0.03, V 3 e = -0.03. 

We note that the volatility surface induced by the T t — t regime is able to produce a negative 
at-the-money (ATM) skew, which is typical of equity call options. However, without a stochastic 
time-change, the implied volatility surface will not exhibit a smile [14] . 

5.2 TC-FMR-SV: Levy Subordinates 

Next, we consider a regime where the random time-change T 4 is given by a Levy subordinator. 
The jumps of our prototype Levy subordinator will be modeled as a compound Poisson process. 
Specifically, we consider 

T/:= 7 t + ^&, (48) 
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where 7 is the drift of the Levy subordinator, N" is a homogeneous Poisson process with jump- 
arrival intensity a and the £j are i.i.d. random variables with exponential distribution ~ £(v) 
and mean E[£j] = l/rj. As noted in section l2.3.1l the Levy measure v{ds) of a compound Poisson 
process can be written as the product of the (net) jump-arrival intensity a and the distribution 
F^(s) of the i.i.d. jumps. In this case 

v{ds) = aF s (ds), F 6 (s) = 1 - e~ r ' s . (49) 

Using equations (J7J and (|49|) . we calculate the Levy exponent 0(A) of a as 

0(A) = 7 A + -^-, (A >-??). 

A + 77 



For a Levy subordinator expectations (|40)) and (|41l) reduce to 



E, 



where, for a compound Poisson process with exponentially distributed jumps, </>L°^ and are 
given by 

aAl 0) , m . m oAJPaP aA^ 



j«J) = ~a<°> + — <£ (1) = tA^ - ^ " + — (50) 



The coefficients Ci 0) (i), given by (|4"T)) . are unaffected by the choice of random time-change. Hence, 
the approximate price of a European call option P^ (t, x) + y/e (t, x) can now be calculated 
using d37]) and (|35|) . 

Figure [2] plots implied volatilities I versus LMMR in the TC-FMR-SV regime in which T t x is 
given by (|48[) . The parameters used in figure [2] are 

r = 0.00,o- = 0.34, V£ = 0.03, V3 = -0.03, 77 = 0.10, a = 0.75, 7 = 0.25 

We note that the implied volatility surface of figure [2] exhibits an ATM skew as well as a true smile 
with implied volatilities rising at the largest strikes. The strong skew and smile are particularly 
noticeable at shorter maturities. This is consistent with the findings of [TS], where it was noticed 
that a model for the underlying asset St must contain jumps in order for the induced implied 
volatility surface to capture the steep skew and strong smile of the empirically observed implied 
volatility surface for short-maturity options. 



5.3 TC-FMR-SV: Absolutely Continuous Time-Change 

Recall that an absolutely continuous time-change T t 2 is of the form (JSJ). As an example, we 
consider Z t to be the classic Cox-Ingersoll-Ross (CIR) process and the rate function to be the 
identity V(z) — z. We have 

dZ t = k(0 - Z t )dt + Y,\/z~ t dW t z , 

T 2 = [ Z s ds, (51) 
Jo 

where is a Brownian motion under P. Here k > is the rate of mean-reversion of the CIR 
process and 8 > is the long-run mean. We shall refer to E > as the "vol of vol" since E 
controls the volatility of Z t , which in turn contributes to the volatility of X T 2. We shall enforce 
the condition 2k® > E 2 so that the CIR process Z t remains strictly positive for all time (see [24], 
Chapter 6). 
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In order to compute option prices in the absolutely continuous time-change regime, we need 
to know the Laplace transform of T 2 . This is a classical calculation, which can be found in [53]. 
Here, we simply state the result 



L(t, z, A) = E z 



-AT? 



-K,eu(t)-zv{t) 



. -K 2 \ 



U(t) 
V(t) 



E 2 



loe 



2<ye 



( 7 +«)t/2 



(7- k) + eT*(7 + k) 



2A (eT* - 1) 



(7-K) + eT*(7 + K)' 



7 = Vk 2 + 2X 2 A. 

For an absolutely continuous time-change T 2 , expectations (|40|) and (|4Tj) reduce to 



E, 



_ A (0) T 2 



(m,a£°)), i,[(-A«T t 2 ) 



A (0) T 2 



= 9, 



^(M,A«a)[ =ALO)/A£i) (52) 



The above expectations, along with expression (|4"T)) for cd, ^), are enough to calculate the ap- 
proximate price of a call option pW(i, x, z) + ^/ePW(t, x, z) using (|37|) and (f38l) . 

Figure [3] plots implied volatilities / versus LMMR in the TC-FMR-SV regime in which T 2 is 
given by (|5T|) . The parameters used in figure [3] are 

r = 0.00, a = 0.34, V£ = 0.03, = -0.03, k = 1.00, 6 = 1.00, S 2 = 2.00, z = 2.00. 

We observe that the implied volatility surface in figure [3] exhibits an ATM skew as well as a slight 
smile effect. Though, neither the skew nor smile in figure [3] is as pronounced as in figure [2] where 
the stochastic time-change is given by a Levy subordinator T t r. 

5.4 TC-FMR-SV: Composite Time-Change 

Finally, we consider a composite time-change Tf = T^ 2 where T} is the Levy subordinator de- 
scribed by equation flU} and T 2 is the absolutely continuous time-change described by equation 
([5Tj) . In this regime expectations (|40|) and (|41"T) can be found by replacing A^ "* and A^ 1 ^ in (|52|) 



by (pS and 4>l} from equation ([J 

In figure U we plot implied volatility / induced by the composite time-change T t 3 as a function 
of LMMR. The parameters used in figure [4] are 

r = 0.00, a = O.M,V 2 e = 0.03, Vi = -0.03, 7 = 0.05, a = 0.50, r] = 0.50, 
K = 2.00, 9 = 1.00, S 2 = 4.00, z = 4.00. 

Once again, we observe an ATM skew and strong smile at the shortest maturity, with these features 
diminishing for longer maturities. 



6 Summary and Conclusions 

In this paper we introduce a class of TC-FMR-SV models. The key features of our modeling 
framework are: 

1. We are able to include jumps in the price process of the underlying asset. 

2. We can incorporate multiple factors of stochastic volatility, which run on different time 
scales. 
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3. We are able to account for the empirically observed negative correlation between asset returns 
and volatility (the leverage effect). 

Some of the main results of our analysis are: 

1. We provide simple formulas to calculate the approximate price of any European option. 

2. By combining different time-changes, we are able to produce a wide array of implied volatility 
surfaces. 

Overall, we feel that the flexibility provided by the TC-FMR-SV framework and the analytic 
tractability it provides, merit continued research in this area. A logical next step, for example, 
would be to incorporate default of the underlying asset into our class of models, as done in 
[2T] , Additionally, characterization of the implied volatility surface through an expansion I e sw 
J<°> + v^ (1) would be useful. 

Thanks 

The authors of this paper would like to thank Jean-Pierre Fouque and two anonymous referees 
for their thoughtful comments on this work. Their suggestions have greatly improved both the 
quality and readability of this paper. 



A Poisson Equations and the Fredholm Alternative 

The purpose of this appendix is to explain why centering condition (1261) is necessary in order for 
the Poisson equation (|25|) to admit a solution. To begin, we consider an ergodic Markov diffusion 
Y t that lives on R, has invariant distribution Fy(dy) = p(y) dy and whose infinitesimal generator 
and adjoint are given by 



p{y)d y 



-dyH{y)+ d yy^— 



From the Kolmogorov forward equation, the density p(y) satisfies £ Y P — and is given by 



p(y) = c 



<? 2 {y) 



exp 



MA 



dz 



where C is a constant such that J p(y) dy = 1. 

We consider the following Poisson problem: find f e C 2 (K) such that 







(53) 



km \^(y)\<\y\P 

y— >±oo 



for some real p < oo. 



Multiplying (£y\f , (?/)) by p(y) and integrating with respect to y we find 



p{C Y ^) dy = Cexp 

t 



p(y)d v ^(y) 



* (CyP) dy 



OO 



where we have used integration by parts and CyP = 0. Hence, 



(g) ■= 



^ 2 (y) 



pgdy = - 



p(y)dy*(y) 



(54) 



In section [3~2l we considered Poisson equations with respect to the operator CS 2 ) = £ y , the 
infinitesimal generator of the volatility-driving process under the physical measure P, which 
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we chose to be an OU process. Under the physical measure the OU process YJ 1 has an invariant 
distribution Fy ~ N(m,v 2 ). In this setting, p(y) in equation (fMf asymptotically behaves like 
~ e~ y . Thus, if we restrict ourselves to solutions ^(y) of (|53|) that have at most polynomial 
growth as y — > ±00 then the right-hand side of is zero. Hence, a necessary condition for the 
solvability of (|53l) becomes (g) — 0, which is precisely the centering condition given in (f2l>)) . 

We have established that (g) = is a necessary condition for the solvability of (|5^|) . It turns 
out (g) = is also a sufficient condition for ()53[) to have a solution. The Fredholm alternative 
states that one of the following is true: 

1. Cy^ + .g = has a unique solution (i.e. Cy is invertible) or 

2. £y*i> = has a non-trivial solution, in which case Cy^+g = has a solution if g _L Ker (Cy). 

For the OU process with infinitesimal generator £y, we note that C Y ^ = has a non-trivial 
solution - namely ^(y) = 1. Hence by the Fredholm alternative £y , I' + .9 = has a solution if 
g _L Ker (jCy). Since we have Ker tC Y ) = {p(y)}, the statement g _L Ker CCy) is equivalent to 
the centering condition (g) — 0. The following (formal) solution can easily be checked 



We refer the reader to section 6.6.3 of [12] for a detailed exposition on Poisson equations and the 
Fredholm alternative. 
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Figure 3: Implied volatility surface induced by an absolutely continuous time-change. 
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